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, We explore the renormahzation group flow of quartic perturbations in the low- 
CN ■ 

^ . enegy theory of graphene, in the strong Coulomb coupling and large- limits, where 

D ■ is the number of fermion flavors. We compute the anomalous dimensions of the 
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quartic couplings u up to leading order in 1 /N and find both relevant and irrelevant 
directions in the space of quartic couplings. We discuss possible phase diagrams and 
relevance for the physics of graphene. 
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^ : I. INTRODUCTION 

I Graphene, a single layer of graphite, has lately attracted increasing attention, especially 



ipn] 

Q I since its experimental realization This material presents a number of unusual electronic 
^ ' features, such as a Landau level structure that gives rise to half-integer quantum Hall effect. 
^ , In general, those properties can be traced back to the low-energy spectrum, governed by 

[ two-component massless fermions that come in four flavors: two due to electronic spin 

and two from degenerate Dirac points in the band structure Q]. The quasi-relativistic 
^ ' electronic spectrum is characterized by a velocity v ~ c/300 <^ c, where c is the speed of 
light. However, the Coulomb interaction, which is instantaneous for all practical purposes, 
breaks this emergent Lorentz invariance. In particular, the electron velocity becomes scale- 
O ■ dependent j^. 

It was noted recently 0] that a generalizaton of the low-energy theory that describes 
^ , graphene possesses a quantum critical point. Reference Q considered a model of N species of 
2 + 1 dimensional two-component massless Dirac fermions interacting through a 3D instan- 
taneous Coulomb interaction with a coupling constant g. It was then shown, by analyzing 
the renormalization group (RG) flow of the fermion velocity, that for sufficiently large A^, the 
strongly coupled limit g'^N/v — oo is a quantum critical point, characterized by a dynamic 
critical exponent z = 1—8/ {n'^N) + 0(1/A^^). It was also argued that real graphene is close 
to this critical point for a large momentum window. 

In the present work we further explore the quantum critical properties of this theory by 
studying the RG flow of various quartic perturbations of the Thirring ^ and Gross-Neveu 
[gI varieties. These four-fermion interactions are irrelevant at weak coupling, but obtain 
nontrivial anomalous dimensions of order 1/A^ at the strong coupling flxed point. If a four- 
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fermi coupling becomes relevant below some A^, one may expect dynamical gap generation, 
in which case the system becomes an insulator, or that the system would flow to another 
fixed point. 

The question of whether the Coulomb interaction makes graphene an insulator has been 
considered previously. In Ref. Q it was found from solving the gap equation with a screened 
Coulomb interaction that, for N < 8/n ^ 2.55, a gap opens at sufficiently large coupling. 
For N > 8/n the system remain gapless at all couplings. This was confirmed in Ref. [s] using a 
similar approach. However, as the approximations employed in these works are uncontrolled, 
one would like to explore alternative approaches to this problem. The RG analysis in this 
paper is one of them. (Instability toward ferromagnetism has also been considered (qI.) 

We implement a Wilson-Fisher RG transformation and find the anomalous dimensions of 
the various couplings. Throughout the paper the Coulomb parameter A = g'^N /32v is kept 
finite and the limit A — oo is discussed at the end. In this limit we identify both relevant 
and irrelevant directions in the parameter space of the quartic couplings and discuss possible 
scenarios, including gap generation and the the fiow to a new non-gaussian fixed point. At 
the end of our results section we address the relevance of our findings for real graphene. 

The paper is organized as follows: the model and Feynman rules are presented in section 
II, the details of the calculation are outlined in section III, the results are discussed in section 
IV and the conclusions presented in section V. 



II. MODEL AND FEYNMAN RULES 



Consider a model of A^/2 flavors of four-component relativistic fermions (or flavors 
of two-component fermions) ipa interacting through an instantaneous Coulomb interaction 
and through generic U(A^/2)-symmetric quartic interactions. The corresponding Euclidean 
action is 

Se = - j dt(fx {iJal%i^a + viJal%i^a + ^^Q^aTVa) + ^ / dtd^x{d,AoY 

2A^ 



where the 7^, /i = 0, 1, 2 (we shall also use Latin indices for the spatial directions) are Dirac 
matrices satisfying a Euclidean Clifford algebra {7^,7^} = 26^''', and 7^ = 'j^'j^'y'^. We can 
choose, for example, a representation of this algebra in which 



where a* are the Pauli matrices. In the real world A^/2 = 2, corresponding to two spin polar- 
izations of the electrons. The two two-component spinors near the two valleys are combined 
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into a four-component spinor. The action (12.11) is invariant under spin rotations, but is not 
invariant under rotations in the valley space. ^ Notice that the four-fermi terms are rotation- 
ally invariant, but not Lorentz invariant: the latter is broken by the Coulomb interaction. 
Apart from the four-fermi terms written in (12.11) . one can also introduce an independent set 
equal to 7^ times the vertices already included, where 7^ is linearly independent from the 
other 7-matrices and {7^,7^} = 0. One possible choice is 

^-(:;)- 

Such matrices, together with the ones already considered, form a complete basis for the 
algebra of 4 x 4 Dirac matrices. We restrict for the moment to the action (12.11) . but we shall 
comment on these extra vertices in our discussion section. 

We shall use the large- limit to do calculations at large values of the coupling g, so 
will remain arbitary (but large) until the end. Notice that, unlike the fermionic degrees 
of freedom, the Coulomb field Aq lives in 3 + 1 dimensions, which is reflected in its kinetic 
term. In the strong coupling g ^ 00 limit this term disappears and the Aq propagator is 
dominated by quantum corrections coming from the fermion loops, as we shall see below. 



Go(p) 
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FIG. 1: Feynman rules (see text for details). 

In 2 + 1 dimensions the naive dimensions of the fields and couplings are as follows: 
[ip] = m, [Aq] = m, [g] = m°, [uj] = m~^. 

The Feynman rules are as in Fig. [T], where the fermion propagator is 

G,{p) = ^, (2.4) 

where p = {pojP), p being a 2D momentum vector. Here and in the rest of the paper, we 
use the notation ^ = ■y^po + ■ p and p"^ = Po + v'^\p\'^- 



The two degenerate Dirac points in the band structure of graphene are sometimes referred to as vaheys. 
Thus, rotations in valley space correspond to operations that mix the upper and lower components of the 
four-component Dirac fermions of our theory. 
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The boson-fermion interaction vertex is 270, and the quartic vertices are —^Tj^Tj, with 

Tj E {l,7o7i,7o75,7o,7i,75}- 

Finally, the bare boson (Aq) propagator is 

DM = ^ (2.5) 

In the large A^, finite g'^N limit we must resum the fermion loop contributions to this 
propagator (see Fig. [2]). 

D{p) 

= ^WV + ^\/\A K/V + K/\i V\/V + 



FIG. 2: Resummation of the one-loop self-energy contribution to the boson propagator. 



q-p 




FIG. 3: One-loop self-energy contribution to the boson propagator. 



Such resummation (which is equivalent to the random phase approximation) results in 
a dressed boson propagator that we shall use in the rest of this work, as shown in Figs. [2] 
and [31 
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(2.6) 



where A = g^N /{32v) 



III. CALCULATION OF THE RG FLOW 

We implement a Wilsonian RG procedure whereby we integrate out modes in the mo- 
mentum shell Ai < p < Aq and then rescale the coordinates as well as the fields. The RG 
equation has the form 

^ = (U - M(A))u(p) (3.1) 
amp 

where [u]j = Ui and M(A) is a 6 x 6 matrix. The 1 in the right-hand side comes from the 
naive dimension of Uj. Furthermore 



M(A) = M,(A) + M^y(A), 



(3.2) 
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where comes from the vertex corrections and M^j comes from wavefunction renormal- 
ization. 



A. Vertex renormalization 



To the first nontrivial order in finding the RG flow of Uj entails computing the 

diagrams in Fig. H] using the Feynman rules described in the previous section. The loop 




FIG. 4: Contributions to the RG flow of uj coming from vertex renormalization. a and b are flavor 
indices and F is a generic four-Fermi vertex. The dashed line is intended to show the difference 
between the loops in the various diagrams. Exchange of dummy indices leads to the same diagrams 
and results in an overall factor of 2. 



integrals have the same form for all three diagrams, namely 

-1 



16 /"^o d^q qf^q" A ( ^ ^ \q\ 



NJm (2vr)3 q' \q\ 



In 



Ao 
Ai 



X Jo (A) for /i = z/ = (3.3) 



X 



ioi n = u = i 



and vanish for ^ ^ v. The integration region is spherically symmetric. Furthermore, notice 
that the integrand is independent of the azimuthal angle, and that the remaining radial and 
angular integrals are decoupled because |g| = 



A/g^VT^-t^, where t = cos 6 and 9 is the 
polar angle, making the calculation straightforward. We have used the following definitions: 



/o(A) 



dt- 



A 



71 2VA2^^^ 



a + Ta^^ 



A (A) 



dt '^^^2 ^ 



^1 l + AVT^ 
In the large A limit, A ^ cxd. 



A Av^ 



:ln 



A + yP^ 



while 



/o(A) -2 + 21n 



A + v^ 



(3.4) 
(3.5) 

(3.6) 
(3.7) 
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The divergence of /o(A oo) is associated the the fact that the Coulomb interaction is 
unscreened at zero momentum and nonzero frequency. As we shall see, however, all such 
divergences cancel out in the j3 functions for Ui. 

With the above identities at hand it is not very difficult to compute the contributions 
of the diagrams in Fig. H] to the f3 functions. The only significant step that remains to 
be discussed is straightforward Dirac algebra to find the products of a small number of 7 
matrices (coming from the vertices and the fermion propagators). For instance, for r4 = 70, 
the calculation of the diagram in the middle of Fig. H] entails computing 

7o7m7o7i.7o = 7o for = = (3.8) 
= —70 ioT fi = u = i (3.9) 

where we only need the case /x = because the integral (13.31) is otherwise zero. The 
contribution of this diagram to the fiow of M4 is therefore 



AU4 = U4 



In 



Ai 



(/o - h) 



(3.10) 



where we have included the factor of 2 coming from an essentially identical diagram where 
the boson propagator appears connecting the fermion lines on the right instead of on the 
left. For this specific vertex the first and third diagrams in Fig. H] cancel out exactly due 
to a sign difference in one of the fermion lines. In this particular case only diagonal fiow is 
generated, i.e. this vertex does not contribute to the fiow of the others. In general, however, 
there is operator mixing: each vertex contributes to its own fiow but also to other vertices. 
Our result turns out to be 
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(3.11) 



B. Wavefunction renormalization 



Upon integration of the high-momentum degrees of freedom, new coordinates and fields 
are defined, that are rescaled versions of the original ones. Doing this explicitly results in 
the following action: 

5eff = - y dt'd\' {Z^'bf^al'd'oiJa + Z^W^a^'d^^a + •.•} (3.12) 

where the subscript "eff" indicates that this action describes modes with momenta \p\ < 
Ai = 61 Aq. Here t' = h^t x' = bix, < bo, hi < 1, and ... contains all the terms in the action 
that are not quadratic in the fermion fields. We shall renormalize the field as 

^„^< = Zo-nrVa (3.13) 
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(3.14) 



and enforce the nonrenormalization of the fermion velocity by requiring 

Z^\%' _ Zoh _ ^ 

which means that the integration shells should actually be non-spherical and have radii 
related by 

(3.15) 



bo = —h 
^1 



However Z^/Zi = 1 + 0{1/N), therefore the deviation from spherical symmetry is small and 
is not important to the order of 1/A^ we are considering. Thus we can compute the loop 
integral, assuming that the integration region is spherically symmetric. 

D{k) 




«7o 



Go{p-k) *7o 



FIG. 5: One-loop self-energy correction to the fermion propagator. 

The one-loop self-energy correction to the inverse fermion propagator was computed in 
Ref . B. In our notation, it corresponds to 



1 - AZq = 1 



:[/o(A)-/i(A)]ln 



Ao 
Ai 



Zfi = l-AZi = l + -— Jo(A)ln 



Nit' 



Ai 



(3.16) 
(3.17) 



The correction to the fermion propagator is suppressed by 1/A^, which is another reason to 
work in the large limit. Without this small parameter the diagram of Fig. ([5]) would be 
of order one. 

To see how the wavefunction renormalization affects the running of the four-fermi cou- 
pling, we write down a generic quartic term in the action 

2 

(3.18) 



/TV 

u / dtd?x j ipqTipa 

^ \a=l 



Rescaling of the momenta and the fields as discussed above leads to a renormalization of 
the coupling constant u according to 

u^u' = bfb^\u + Au) (Zq^^%^^) ^ = bi{u + Au + u{AZo + AZ^)) + 0{u'^/N) (3.19) 



where we have also included the contribution Au from the vertex renormalization, and where 
AZq + AZi comes from wavefunction renormalization. 
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From Eqs. (13.161) and fl3.17p we find 

M./(A) = 



iVTr" 



■[2/o(A)-/i(A)] 



(3.20) 



IV. RESULTS AND DISCUSSION 



Combining Eqs. (13.111) and (I3.20p . we find 
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(4.1) 



Notice that the terms containing /o(A) cancel out, as anticipated, so the limit A 
finite. The eigenvalues yield the anomalous dimensions: 7^ = 
The two operators with lowest anomalous dimensions, — 40/A^7r^ 



00 IS 

8 [5,-1,-1,-1,-1,5]. 



are 



Including the vertices of the form 'y^Tj ® T^Tj, that were mentioned in the introduction, 
one can repeat the calculation and obtain as a result a copy of the above matrix. This is 
simply a consequence of {7^, 7^} = 0. There is no mixing between these new vertices and 
the ones considered in our calculation. The list of the most relevant operators may then be 
extended to include 



2(V^a7^75^a)^ - (V^a7SV^a) 



(4.3) 



whose the anomalous dimensions are also — 40/A^7r^. 

In the large limit, the anomalous dimensions are small, hence all four-fermi operators 
have dimensions close to 4 and are irrelevant. However, as one decreases the dimensions 
of some operators decrease as well. Naively, at 



N < iV,H 



crit 



40 



4.05 



(4.4) 



the operators in (14. 2 p and (14. 3 p would have dimensions less than 3, and become relevant. Of 
course, this is only an extrapolation of our leading-order result to finite A^. Nevertheless, it 
is quite interesting that A'^crit is close to the real-world value of A^ = 4. (For comparison, for 
the Lorentz-invariant Thirring model the critical number of four-component Dirac fermions 
is quoted as 6.6(1) [l3], to be compared with A'crit/2 ~ 2.) Due to the limitation of our 
calculations, we cannot determine whether the exact A"crit is smaller or larger than 4. 

If A'crit < 4, then the physics at A^ = 4 and infinite coupling is governed by the strong- 
coupling fixed point discussed in Ref. Q. If A'crit > 4, then there are two further possibilities 
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for = 4, A = oo. It may turn out that the system develops a gap, and becomes an insulator. 
In this case one expects a bifermion operator to have a nonzero expectation value, breaking 
a discrete symmetry. Unfortunately, simply from the form of the operators (14.21) and (14. 3 p 
it is not possible to conclude, with definiteness, which of the discrete symmetries will be 
spontaneously broken. Another possibility is that the coupling flows into a new stable fixed 



point (a similar situation was discussed in Ref. In this case the system remains gapless, 
but with a new dynamic critical behavior. This case is illustrated in Fig. [6l where we also 
show the P function of the coupling u. 




FIG. 6: Phase diagram and beta functions (Inset), in one particular scenario. 

In real graphene, on the other hand, A is large but not infinite, so in the above discussion 
we should take into account that the anomalous dimensions are functions of both and A. 
In real graphene one has 

32eohv ^ ^ 

where = 4 and v = 10^ m/s, so that the anomalous dimension of the operators (14. 2 p 
and (14.31) is 7 ~ —0.72. Thus, all four-fermi operators are irrelevant, at least to the order of 
we are considering. It is conceivable, however, that higher-order corrections will push 
7 to be below —1. On the other hand, for graphene on a Si02 substrate with dielelectric 
constant e = 5.5, the coupling A is reduced by a factor of 2/(1 + e), and 7 ~ —0.45, 
substantially above —1. In this case, one can conclude that the four-fermi interactions 
are irrelevant (and will become more irrelevant as the fermion velocity v increases in the 
infrared) . 
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V. CONCLUSION 

We have studied the RG flow of various four- Fermi couphngs u in the low-energy theory of 
graphene, generalized to include different fermion flavors. We computed the anomalous 
dimensions of the various couplings to the first nontrivial order in In the limit of 

infinitely strong Coulomb interaction, the operators with lowest dimensions become relevant 
for < A^^crit, where iVcrit is estimated to be 4. In real graphene with = 4 but at finite 
Coulomb coupling, the four-fermi interactions are irrelevant, at least to the leading nontrivial 
order in 1/A^. 

It would be interesting to further investigate the phase diagram of our model. One can 
try to push the calculations to another order in 1/A^. In addition, one should try to perform 
numerical simulations of the model. The most interesting values of A^ where nontrivial 
phases may exist, as indicated by our calculations, are lower values like N = 2 and A^ = 4. 
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Note added: After this work was completed, the authors learned of Ref. 
the running of four-fermi interactions is also considered. Our lowest anomalous dimension 
coincides with the value found in Ref. 12, but our values of other anomalous dimensions 



disagree with Ref. We thank Igor Aleiner for bringing Ref. 
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